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Using torrefied biomass, or biocoal, as a solid combustion fuel provides opportunity to introduce a sus¬ 
tainable feedstock into the energy market. The decomposition of biomass to biocoal is a very complex 
chemical process, yet detailed understanding of the mechanistic changes that occur can allow for sys¬ 
tem control and optimization. Torrefaction of lignocellulosic biomass has been studied extensively, and 
kinetic models have been developed to attempt to describe the observed behavior during Thermal Gravi¬ 
metric Analysis (TGA) experiments, using very slow heating rates. The goal of this paper was to develop 
a baseline semi-empirical model for torrefaction of aspen wood using very fast heating rates (1000 °C/s) 
and GC-MS analysis. This fast heating rate was adopted in order to uncouple heat transfer and chemical 
kinetics. 

Having direct information about the specific chemical evolution during torrefaction affords mecha¬ 
nistic insights into torrefaction, along with opportunities for improved process design and control. The 
proposed kinetics model involves three sequential reactions in which initially the parent hemicellulose 
material produces gaseous products (CO and C0 2 ), water, and organic acids (formic and acetic) and a 
solid intermediate. The next reaction further degrades the solid intermediate into the stated products. 
The final reaction further generates the products along with the final bio-coal product. We found the 
sequential reactions to have characteristic times of 0.99 min, 2.14 min and 39.45 min respectively at 
300 °C. Improvements to the models and routes for future investigation are also discussed. 

© 2013 Elsevier B.V. All rights reserved. 


1. Introduction 

Coal-fired utility boilers are the most prominent electricity gen¬ 
eration systems, producing 40-50% of consumed electricity across 
the globe (42% in the U.S.) [1,2], Coal exists in great abundance 
around the world, and is a valuable solid fuel because of its high 
energy density and low cost. Coal does, however, have well-know 
environmental issues such as the production of CO, CO2, SO x , 
NO x , particulates, mercury, and poly-aromatic hydrocarbons to 
name a few [3-5], In the current age of increasing environmental 
awareness, these issues inherent to coal combustion must be con¬ 
tinuously addressed and controlled through the use of additional 
high-cost process operations. 

A possible solution to environmental concerns is replace¬ 
ment of the coal by a different feedstock that is environment 
conscious-sustainable biomass. Because the coal-power industry 
is so prominent, it is important, if not necessary, to use a drop-in 
replacement in existing infrastructure. Initially, raw regenerative 
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biomass was investigated as a possible solution, and demonstrated 
for commercial purposes [6,7], The combustion of biomass inher¬ 
ently has a preferable emissions profile because it produces no net 
C0 2 emission due to a biogenic carbon cycle, and does not contain 
sulfur and mercury [8-11], However, some issues exist with raw 
biomass processing and combustion, including high bulk volumes, 
moisture content, low calorific value, material hydrophobicity, pro¬ 
duction of undesirable tars, and high energy and time requirements 
for size reduction [12-14], One solution to address these issues 
and still utilize biomass based fuel, has been the implementation 
of torrefaction as a pretreatment to produce biocoal. 

Biomass torrefaction is a method of mild thermal treatment 
in the temperature range of 200-300 °C under inert atmosphere 
that, through a complex coupled chemical-kinetic-heat-and-mass- 
transfer process, converts the biomass into biocoal. In this process 
mostly the hemicellulose fraction of biomass degrades, with small 
contributions from the lignin and cellulose fractions [15-17], The 
products of torrefaction are 70-90% solids (ash and biocoal), 6-35% 
liquid (mostly water with some fatty acids, phenols, acids, alcohols, 
furans, ketones) and 1-10% gas (CO, C0 2 , H 2 ) [16,18], 

There are major advantages to biocoal when compared to raw 
biomass, one of the most prominent being reduction in energy and 
time requirements for particle size reduction. The energy required 
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Table 1 

Ultimate and proximate analysis of coal, biocoal, and various raw biomass species [36-40], 


Ultimate analysis Proximate analysis 

C(%) H (%) 0(%) N (%) S(%) FC(%) VM(%) Ash (%) HHV (MJ/kg) 


Corn Stover [36] 42.6 

Softwood [37] 52.1 

Hardwood [37] 47.8 

Switchgrass [36] 42.0 

Lignite [37] 63.9 

Sub-bit. Coal [38] 70.0 

Biocoal [39,40] 70.7 


FC: Fixed Carbon, VM: Volatile Matter, 


5.1 36.5 0.8 

6.1 41.0 0.2 

5.8 45.8 0.1 

5.2 34.0 0.7 

5.0 24.5 0.6 

5.2 23.6 1.0 

5.6 20.1 0.6 


HHV: Higher Heating Value 


15.4 

28.1 

11.3 

12.9 


32.3 


69.7 

70.0 

88.2 


49.5 
48.0 

63.5 


16.3 

20.1 


16.3 

17.6 

25.9 

28.5 


to pulverize biocoal is approximately one tenth that of biomass [ 1 2 ] . 
This is due to the decomposition of hemicellulose which changes 
the viscoelastic properties of the feedstock [19], In addition to this 
feedstock improvement, torrefaction also increases the calorific 
value of biomass up to 40% by removing bound feedstock mois¬ 
ture and reducing the mass by up to 30-50% [20,21], Finally, the 
oxygen to carbon ratio is improved by significantly reducing the 
oxygen content, while the carbon fraction is increased accordingly 
[21], Mostly importantly, torrefied biomass (biocoal) has similar 
physico-chemical properties to traditional coal and can be con¬ 
sidered a premium drop-in replacement, requiring no addition 
upgrades to existing utility boilers. Table 1 compares the properties 
of coal, biocoal, and various raw biomass species. 

Regardless of the feedstock or energy production method, how¬ 
ever, a key aspect of electricity-generation systems is their high 
feedstock throughput [7,22], As demands for energy are projected 
to increase consistently by approximately 50 quadrillion BTUs 
every 5 years starting from a measurement of 505 quadrillion BTUs 
in 2008 [23], a detailed understanding of the entire energy gen¬ 
eration process is required for process optimization and design. 
Because of similarities between the properties of biocoal and tra¬ 
ditional coal, most processing steps are similar to those already 
in practice with exception of feedstock production. For this rea¬ 
son among others, some emphasis has been placed into developing 
models and detailed understanding of the production of biocoal, or 
more specifically on torrefaction technology. Recent review arti¬ 
cles have summarized several existing proposed wood torrefaction 
models based on Thermo Gravimetric Analysis (TGA) including a 


single step global model [24], a three parallel reaction Shafizadeh 
and Chin model [25], a two step consecutive Di Blasi-Lanzetta 
model [24,26,27], and a two parallel Broido-Shafizadeh model [24], 
Despite these efforts there are many competing ideas about the 
complex chemical pathways, with no definite solution [17], 

The present work focuses on proposing a mechanistic reaction 
scheme with a semi-empirical model, using different techniques 
than traditionally employed. The proposed models were developed 
with very high sample heating rates (1000°C/s), and direct Gas 
Chromatography/Mass Spectroscopy (GC/MS) analysis. Very high 
heating rates during model development are critical in order to 
uncouple heat transfer from chemical reaction and thus study the 
true intrinsic reactions present during torrefaction. 

2. Materials and methods 

2.1. Biomass preparation 

Debarked aspen samples were obtained from Dr. Christopher 
Webster in the School of Forest Resources and Environmental Sci¬ 
ence at Michigan Tech. Samples were debarked and then dried from 
their original moisture content at 105 °C in a drying oven until a dif¬ 
ference in mass was no longer detectable. In addition, samples were 
milled to smaller particles using a small scale Thomas Wiley® knife 
mill (NR. 3557524 359264). Following size reduction, a W.S. Tyler 
Rotap (model RX-29, serial 9774) equipped with a range of sieve 
trays was used to obtain samples with known dimensions. Sizes 
larger than 32 Tyler Mesh (500 microns) but smaller than 28 Tyler 


Pyrolysis Gas Chromatography Mass Spectroscopy 



Sample 

Probe 


Fig. 1. Experimental flow diagram. 
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Mesh (599 microns) were considered for this study. The fibers were 
observed to have approximate shapes of cylinders. The wood fibers 
had diameter of approximately 0.5 mm and lengths ranging from 
1.0-3.0 mm. 


2.2. Torrefaction-GC/MS experiment 

Torrefaction of aspen samples in this work was performed at 
300 °C for 90 min with a micropyrolysis reactor (model 5200HP 
Pyroprobe) manufactured by CDS Analytical. A microbalance (Citi¬ 
zen Scales Inc., Model CM5) with a readability of 1 p,g was used to 
record the sample masses. Fig. 1 shows the experimental apparatus. 

During the experiment, approximately 0.5 mg biomass samples 
were contained within a quartz reaction vial with an inert ultra high 
pure (99.999%) helium atmosphere. Sample heating (1000 °C/s) is 
achieved via radiative, convective, and conductive heat transport 
from a platinum filament surrounding the transparent reaction 
vial. Once sample heating begins, the flow of helium gas over the 
sample is diverted to the inlet port on a Trace GC Ultra (Model 
K8880181) gas chromatograph by Thermo-Finnigan (now The- 
moFisher) equipped with a 30 m RXI-5MS fused silica (low polarity 
phase, Crossbond® 5% diphenyl/95% dimethyl polysiloxane) capil¬ 
lary column by Restek. The GC column was maintained at 300 °C 
during experimental trials, and had a constant 1.5 ml/min (superfi¬ 
cial velocity of 50.9 cm/s) gas throughout. The GC was connected to 
a quadrupole mass spectrometer (Trace DSQII, Thermo-Finnigan), 
which recorded the ion fragment intensities over time. The MS 
was set to monitor a range of molecular mass from 1-350 m/z. 
The experimental conditions mentioned above provided an envi¬ 
ronment that minimizes interactions of compounds with system 
components (GC column, etc.). 

In addition to the experimental procedures above, great care 
was taken to evaluate complicating factors due to heat and mass 
transfer with the sample, column interaction effects, and flow dis¬ 
persion. These evaluations combined both theoretical calculation 
and using pulse-inputs into the GC under conditions identical to 
the experiments. Input-response results are shown in 

Fig. 2, and directly show how the instrument system interacts 
to create the observed signal. These results are discussed in the 
following sections. 

3. Results and discussion 

3.1. Method validation 

The measurement of torrefaction products is affected by sev¬ 
eral factors. A complicating factor of great importance is limitations 
associated with heat transport. Using solutions for unsteady-state 
heat transfer in solids tabulated in work by Carslaw and Jaeger [28], 
it was found that approximately 90% of the surface temperature 
(boundary value) is observed by the center of a finite cylinder (with 
radius of the wood chips) after only 0.24 s. For these calculations 
the maximum effective radius (300 microns) allowable through 


the sieving step was considered for the characteristic length, along 
with a thermal diffusivity of 1.79 mm 2 /s for biomass [29], Detailed 
results and temperature profiles are shown in the supplemental 
material. Given that the time scale for conducting energy transport 
(0.24 s) is much less than the duration of the experiment (90 min), 
heat transfer limitations are considered to be very minor. The radia¬ 
tive heat transport from the heated filament to the biomass particle 
surface was not considered due to the effectiveness demonstrated 
in the slower of the two major heat transport steps, conduction. 

Similar to the analysis addressing heat transport limitations, 
mass transport was investigated assuming a porous media model 
for the torrefied biomass. Using methods similar to heat transfer 
study, the transport is represented by an effective diffusion coef¬ 
ficient, and a value of 1.50mm 2 /s was assumed [30,31], The mass 
transport becomes analogous to the heat transport, and in the case 
of torrefaction for the stated particle size mass transport limitations 
from diffusion is insignificant. 

Finally, interactions of the products of torrefaction with the 
GC/MS instrumentation (GC column resin) were considered. An 
example of signal convolution due to the GC/MS is shown in Fig. 2 
As indicated above in these tests, a GC auto sampler injects a sam¬ 
ple (1 pi) of acetic acid at time zero; and this is represented by 
a delta function input. It is important to note that the amount of 
acetic acid injected was much greater than that observed over the 
time of torrefaction (90 min)-meaning that any observed interac¬ 
tion effects due to the system are accounted for liberally by this 
method. After approximately 62 s, the detector begins to observe 
a significant amount of the acetic acid ion fragment. Based on the 
superficial velocity, the time the gases/vapors take approximately 
59 s to flow through the column. The output, or measured prod¬ 
uct signal, is dispersed in an almost Gaussian fashion, with a Full 
Width at Half Maximum (FWHM) of approximately 1.4 s. These 
results, when compared to the experiment duration of 90 min, indi¬ 
cate that there is negligible interaction with the chromatography 
system and the modeled compounds beyond the initial production 
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rate (<0.5 min). We conclude from these validation analyses that 
the measured response at the mass selective detector is indicat¬ 
ing the actual rates of production of the species in the torrefaction 
chamber of the micropyrolysis instrument. 

3.2. Model development 

Although there are many identifiable compounds from torrefac¬ 
tion, the model in this work was targeted around major compounds 
produced during torrefaction, namely: CO, CO2, water, and organic 
acids [21], Fig. 3 shows a comparison between the MS relative 
intensities obtained for these compounds, represented by a char¬ 
acteristic ion fragment. The different plots represent the same data 
at varying time scales to highlight their similarities and differences 
- namely that their initial production are very similar, but several 
compounds exhibit a major secondary production peak. The char¬ 
acteristic fragments used were chosen based on their abundance 
in the species’ fragmentation pattern: CO (28 m/z), CO2 (44m/z), 
water (18 m/z), acetic acid (60 m/z), and formic acid (29 m/z). Other 
species commonly associated with the 60 m/z fragment (glyco- 
laldehyde and methyl formate) during thermal degradation were 
not detected, and are thought to originate from the decompo¬ 
sition of cellulose structure, which is not degraded significantly 
during torrefaction [32,33]. The carbonyl (29 m/z) fragment which 
is associated with a carboxyl group, was found to more readily 
form carboxylic acids as opposed to other organic compounds 
[21], Finally, because 29 m/z was not observed as a characteristic 
ion fragment of acetic acid, this carbonyl fragment was associated 
directly to the other carboxylic acid detected, formic acid. 

All compounds shown were normalized by the height of their 
first production peak for comparison (MS relative intensities), as 
some compounds exhibit more than one peak during torrefaction. 
The first production of all the compounds appears to be identical 
(characteristically) for all compounds at early times, though differ¬ 
ent as torrefaction advances. Examining these traces on a log-linear 
plot (Fig. 4) at early times (0.5-2.0 min) can aid in understanding the 
kinetics of production at early times. If relative intensity changes 
in a linear fashion when plotted as in Fig. 4, then the initial produc¬ 
tion is described by an exponential (first-order irreversible) decay 
process. From Fig. 4 if we exclude the first 0.5 min and focus on 


the behavior afterwards it is clear that these compounds are pro¬ 
duced with similar characteristics (overlapping lines indicate the 
similar production) and with first order kinetics - i.e. they are pro¬ 
duced simultaneously in the same first order reaction. From this, 
it is proposed that the first reaction in the torrefaction mechanism 
is the decomposition of the biomass hemicellulose to directly pro¬ 
duce the various products. If this is true, the production rates must 
be at their maximum at time zero. This discrepancy with what is 
observed is caused by the small convolution imparted to the signal 
by heat-mass transfer as well as chromatographic influences previ¬ 
ously discussed. It is proposed, then, that the first order exponential 
function that describes times 0.5-2.0 min should be extrapolated to 
describe the first reaction from time zero. 

However, upon examining the production beyond these early 
times, it is noted that all compounds are continually produced yet 
with something other than a first order decay. Namely, the major¬ 
ity of acetic acid is produced beyond this first mechanism. The 
only way for this to occur is through the use of solid intermediates 
between the parent biomass and the final bio-coal. These interme¬ 
diates would then further degrade to produce the observed species 
in subsequent irreversible first-order reactions. It is proposed that 
there are two additional reactions beyond the first, where solid 
intermediates produced from the previous reaction further degrade 
to produce the observed species. 

3.3. Kinetic model and interpretation 

The proposed set of reactions begins with the degradation of the 
parent hemicellulose material (H) into a solid reaction intermediate 
(1?!) and products from reaction 1 (Pj ,) with reaction rate constant 
k\. The second reaction consists of the solid reaction intermediate 
further degrading through a second reaction rate constant (/C2) to 
further produce products (P, 2 ) and intermediate (R 2 ). The third and 
final reaction produces products (P,- 3 ) and the final bio-coal product 
(C) from the second intermediate with rate constant (k 3 ). These 
reactions are summarized in Eqs. (1 )—(3). 


H XR, + ^iPi 

(1) 

Ri^R 2 + ^2fiiPi 

(2) 

Ri-^C+Y^YiPi 

(3) 


where a„ fo, y,- are stoichiometric coefficients for production of 
species i in reactions (1)—(3), and k\ >k 2 > k 3 . Assuming first order 
reactions, the time dependent solution for the rate of production of 
product P, is shown in Eqs. (4)—(6), and the details of this solution 
are shown in the supporting materials. 

^ = a,k,H + + Vik 3 R 2 (4) 

Contribution from 
(reactionl) 

(5) 

(reaction2) 


=a i 0[exp(-k 1 t)] 

+/J, Oipo) [exp(-fc, t) - exp(-fc 2 t)[ 


+YiO(pco [(k2 -fc 3 )exp(-k]t) + (k 3 -ki)exp(-k2t)-(k2 - ki)exp(-k 3 t)[ (reaction3) 


9=k x H 0 , 


k 2 k 3 

v k 2 - kl’ (k 3 -k r )(k 2 - k 3 ) 


(6) 


The reaction variables fci, k 2 , k 3 , ft, y,- and (H) 0 were fit simul¬ 
taneously to all data by minimizing the squared error between 
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Fig. 5. Experimental rate data with model predictions. 


the measured experimental data and the model predictions. The 
objective function to fit the variables is shown in Eq. (8). 


min EE( p y p y) 2 

(8) 

j i=0 


Pij = P(ti,{ki,k 2 ,k 3 ,p j ,y j , Wo }) 

(9) 


where py is the measured data for species j at time i, Py is the 
calculated model value for species j at time i. The reaction rate 
coefficients were found to be k\ = 1.0079 min -1 , k 2 = 0.4664 min -1 , 
and k 3 =0.0253 min -1 . The characteristic times, or inverse reac¬ 
tion rates, were 0.99 min, 2.14 min and 39.45 min respectively. The 
experimental data with the reaction scheme are shown in Fig. 5. 
This figure shows very close agreement between the proposed 
model (labeled “Model’’ in the figure) and the experimental data 
(labeled 44 m/z, 60 m/z, etc.), as well as the contributions from the 
three reactions (labeled “Rxn 1”, “Rxn 2”, “Rxn 3”). 

Each compound shares an identical production according to 
the first degradation reaction, which arises from the normaliza¬ 
tion of the data. During this method, a ,• is forced to a value found 
with the extrapolation mention previously, while the other stoi¬ 
chiometric factors were calculated as relative values. The values of 
these coefficients are shown in Table 2. The model predicts that the 


maximum amount of hemicellulose occurs at time zero, while the 
first intermediate peaks around 1.5 min, and the second intermedi¬ 
ate is maximized after approximately 8 min. It is clearly shown that 
the second reaction has some contributing effects, but the majority 
of sustained chemical production occurs during reaction 3. Per¬ 
haps the solid hemicellulose intermediate at this stage is unstable, 
or readily fragments to a more suitable derivative. For all species 
studied here, the model predictions are in good agreement with the 
laboratory data. This lends evidence to our proposed model, sug¬ 
gesting that the products are produced simultaneously is a series 
of sequential reactions. This idea is shared by other works which 
attempt to describe torrefaction. Table 3 displays a summary of 
some established torrefaction models [17], 


Table 2 

Stoicheometeric factors for each modeled compound. 


CO 

C0 2 

H 2 0 

CH3COOH 

HCOOH 


1.28 

1.28 

1.28 


24.2 

8.65 


60.0 

15.0 
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Table 3 

Torrefaction kinetic models, adapted from [17,24,26,34,41 ]. 


Kinetic model 


Mechanism scheme 



Source 

Three consecutive reactions 


Hpmirpllulnsp 

ki 

Intermediate 1 + Volatile 

This work 



k 2 

- Intermediate 2 + Volatile 




k 3 

Char + Volatile 


One step global 


Hemicellulose 

ki 

y Char + Volatile 

[17,24] 


ki 

y Volatile 

Two consecutive parallel reactions 


Hemicellulose 

k 2 

Intermediate 

[26,34] 




k 3 

■j Volatile 




Intermediate 

k4 

Char 





ki 

j Tars 


Two parallel reactions 


Hemicellulose 

k 2 

Char + Volatile 

[17,24] 




ki 

Volatile 


Three parallel reactions 


Hemicellulose 



[17] 





"i. 

Char 





k l, 2 ,3 

Volatilei 23 


Five parallel reactions (3 decomposition h 

■ 2 oxidation) 

Hemicellulose 

k4,5 

^ Volatile 4 5 

[41] 


When compared to existing models, the work here has sev¬ 
eral novel aspects and insights: (1) moves beyond lumped models, 
to the identification of specific chemical species transients’ dur¬ 
ing torrefaction, (2) the model was developed around observed 
species, and thus direct changes to the biomass instead of weight 
loss curves, (3) very fast heating rates, which decoupled transport 
and kinetic effects arising from temperature ramping. When com¬ 
paring to existing models, perhaps the closest match is the two-step 
consecutive Di Blasi-Lanzetta model. This model has been used 
extensively in other studies in attempt to explain trends obtained 
from TGA studies [26,34,35], In this model, the parent biomass (A) 
forms volatile product group 1 (Vi) and a reaction intermediate 
(B) through reaction rates k v] and kj respectively. After this first 
reaction, the intermediate is further broken into a second group of 
volatile products ( V 2 ) and a solid product (C). These reactions are 
summarized in Eqs. (10)—(11) [26], 

A-Hvi,A-Hb (10) 

bHv 2 , b\c (11) 

4. Conclusions and recommendations 

The new model presented in this work describes the evolution 
of gases (CO and CO2), water, and organic acids (formic and acetic) 
from the torrefaction of aspen. The proposed models were formed 


around the observations of fast product evolution followed by vary¬ 
ing amounts of sustained production at later times. The model 
consists of three sequential reactions with characteristic times of 
0.99 min, 2.14 min and 39.45 min respectively. The data for the 
species measured were well fit by the model. Transient characteri¬ 
zation of the remaining biomass solids could provide a well-defined 
path for model improvements. Techniques such as Fourier trans¬ 
form infrared spectroscopy (FTIR) could be employed to develop 
a key understanding of change in functional groups. In addition, 
convoluting factors such as heat transfer, flow dispersion, and GC 
column interaction should be further addressed or incorporated in 
the model to further refine the key early-time assumptions. 

Appendix A. Supplementary data 

Supplementary data associated with this article can be found, in 
the online version, at http://dx.doi.Org/10.1016/j.jaap.2013.08.010. 
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